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We present a Monte Carlo simulation study of measuring the rare leptonic decay B^ — * (J, + fi~ 
with the CMS experiment at the LHC. The study is based on a full detector simulation for signal 
and background events. We discuss the high-level trigger algorithm and the offline event selection. 
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1. Introduction 

In the standard model (SM) , the decay B® — > 
has a highly suppressed rate 1 of B = 
(3.42 ± 0.54) x 10~ 9 since it involves & b ^ s 
transition and requires an internal quark an- 
nihilation which further suppresses the decay 
relative to the electroweak 'penguin' b — ► 
decay. In addition, the decays are helicity 
suppressed by factors of mf . To date these 
decays have not been observed; upper lim- 
its on these branching fractions are a topic 
of frequent updates at the B-factories 2,3 (for 
the decay 5°, — > /x + /i~) and the Tevatron. 4 ' 5 
Currently the best upper limit is from the 
CDF collaboration 5 with B(B° -> m + aO < 
1.0 x 10~ 7 at 95% confidence limit. 

Since these processes are highly sup- 
pressed in the SM, they are potentially sen- 
sitive probes of physics beyond the SM (see 
Fig. 1). In the minimal supersymmetric ex- 
tension of the SM (MSSM) the branching 
fraction for these decays can be substan- 
tially enhanced, especially at large tan/3. 6 
For MSSM with modified minimal flavor vi- 
olation at large tan /?, the branching frac- 
tion can be increased by up to four orders of 
magnitude. 7 B®,* — > decays can also 

be enhanced in specific models containing 
leptoquarks 8 and supersymmetric (SUSY) 
models without R-parity. 9 

There has been some interest 10 ' 11,12 in 
using the decay mode B° s — > n + to 'mea- 
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Fig. 1. Feynman graphs for the decay B^r^ ~* 
illustrating possible new physics contribu- 
tions. 

sure' the key parameter tan j3 of the MSSM 
and to constrain other extensions of the 
SM. Lower bounds on tan (3 can be obtained 
from B(B® — > /i + /i~) with general model- 
independent assumptions. Since tan (3 is also 
constrained from above due to general princi- 
ples, a lower bound is tantamount to a mea- 
surement of tan/3. 

2. The CMS Experiment 

The Compact Muon Solenoid (CMS) detec- 
tor is well suited for the study of leptonic B 
decays. The main components for this anal- 
ysis are the tracker and muon systems. 

The CMS tracker is an all-silicon detec- 
tor, comprises an inner pixel vertex detec- 
tor plus an outer strip track detector, and 
is immersed in a magnetic field of 4 T. The 
pixel detector provides high-precision mea- 
surements of space points close to the inter- 
action region for vertexing and effective pat- 
tern recognition in the high-track multiplic- 
ity environment at the LHC. The pixel de- 
tector is composed of 1440 modules arranged 
in three barrel layers (at a radial distance of 
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r = 4.4, 7.3, 10.2 cm from the beampipe) and 
four endcap disks (at z — ±34.5, ±46.5 cm 
from the interaction point). The barrel de- 
tector comprises 672 modules and 96 half- 
modules, the forward detector is built of 96 
blades with 672 modules. With a pixel size of 
dfj, x d z — 100 /an x 150 /an, a hit resolution 
of 10-20 ^m is achieved. The strip tracker 
allows the precise determination of charged 
particle momenta. It is arranged in the cen- 
tral part as inner (TIB) and outer (TOB) 
barrel, and in the forward regions as inner 
discs (TID) and endcaps (TEC). The barrel 
part consists of 10 layers (4 in the TIB, 6 in 
the TOB) and 12 layers in the forward part 
(3 in the TID, 9 in the TEC). The pitch in 
the strip tracker varies between 80-180 /jm. 
The material inside the active volume of the 
tracker increases from « 0AX o at pseudo- 
rapidity r\ = to around lXo at \rj\ s=s 1.6, 
before decreasing to w 0.6Xq at \rj\ = 2.5. 

The CMS muon system, incorporated 
into the magnet return yoke, is divided into 
a barrel < 1.2) and forward parts (1.2 < 
\r]\ < 2.4). In the barrel region, where the 
neutron induced background and the muon 
rate is small, drift tube (DT) chambers are 
used. In the two endcaps cathode strip cham- 
bers (CSC) are deployed. In addition, resis- 
tive plate chambers (RPC) are used both in 
the barrel and the endcap region. The RPC 
spatial resolution is coarser than for the DT 
and CSC, but their excellent time resolution 
allows the unambiguous identification of the 
correct bunch crossing. 

3. Event Samples 

Monte Carlo (MC) event samples were gener- 
ated with Pythia 6.227 and passed through 
a full detector simulation based on Geant 4. 
On average five pile-up events were included, 
appropriate for a luminosity of C = 2 x 
10 33 cm-'s" 1 . 

Both signal and background MC event 
samples have been generated as minimum 



bias QCD events. In the signal sample, B s 
mesons decay as B® — > The muons 

are required to have transverse momentum 
p±' i > 3GeV and |<| < 2.4; the B s must 
have p± Bs > 5GcV. 

The background sample contains two 
muons with p± > 3 GeV and \t]\ < 2.4. Their 
separation in azimuth and pscudorapidity 
AR(pn) = \J Acf> 2 + Ai] 2 is required to be 
0.3 < AR(fj,fi) < 1.8. Currently the back- 
ground simulation docs not include muons 
due to hadronic in-flight decays or punch- 
through. It is estimated that this hadronic 
component will increase the background level 
by about 10%. Background events from rare 
Bd, B u , B s , B Cl Ab decays are not included. 

In total 20000 signal events and about 
15000 background events have been ana- 
lyzed. The small size of the background sam- 
ple is the limiting factor of the present study. 

4. Trigger Strategy 

The CMS detector has a twofold trigger 
strategy. The first level trigger, with a la- 
tency of 3.2 /i s, is based on information from 
the calorimeters and the muon system. The 
threshold for inclusive isolated single muons 
is at p±^ > 14GcV, for dimuon events both 
muons must have p± > 3 GeV. The ex- 
pected trigger rates amount to 2.7 kHz and 
0.9 kHz, respectively. The LI trigger output 
rate is at most 100 kHz. 

The high-level trigger (HLT) is a soft- 
ware trigger, running on a large processor 
farm. The HLT reduces the overall trigger 
rate by three orders of magnitude to about 
100 Hz. To fit into the tight time constraints 
imposed by the high input rate, tracking at 
the HLT is sped up by two concepts: (1) 
'Regional seed generation' limits track seed- 
ing to specific regions of interest, e.g., a cone 
around the LI muon candidate direction. (2) 
'Partial tracking' pursues track reconstruc- 
tion only until some criteria are met, e.g., a 
resolution of 2% on the transverse momen- 
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turn. Already with six reconstructed hits, 
both the efficiency and the resolution are 
comparable to the full tracking performance. 

The HLT strategy for B® — > pro- 
ceeds along the following path. (1) Verifi- 
cation of the two LI muon candidates. (2) 
Tracks reconstructed only with the pixel de- 
tector are used to compute a list of possible 
primary vertices, the three most significant 
are retained. (3) Regional track reconstruc- 
tion with up to six hits is performed in cones 
around the LI muon candidates. (4) Recon- 
structed tracks with p± > 4 GeV are paired 
and retained if their invariant mass falls into 
predefined regions for the signal and side- 
bands. (5) The two tracks must have op- 
posite charge and are fit to a common ver- 
tex; the event is retained only when the fit 
X 2 < 20 and the three-dimensional flight 
length Isd > 150 /im. With this selection, the 
event rate was estimated 13 to be < 1.7 Hz. 

5. Offline Analysis 

The offline analysis selection focuses on a sec- 
ondary vertex that is well measured and sep- 
arated from the primary vertex and consis- 
tent with the decay of an isolated B s meson. 

The primary vertex is determined from 
all tracks with p± > 0.9 GeV, because of the 
rather low- multiplicity track environment. 

The two muons must have opposite 
charge, p^ > 4.0, and \rf\ < 2.4. The 
azimuthal and pseudorapidity separation of 
the two muons 0.3 < AR(fj,fi) < 1.2 pro- 
vides a powerful reduction of gluon-gluon fu- 
sion background with both 6-hadrons decay- 
ing semileptonically: The muons of those b- 
hadrons tend to be back-to-back, while the 
signal shows a peaked distribution with a 
maximum at AR(fifi) <~ 1. 

B a candidates are formed by vertexing 
the two muon candidates. The vertex qual- 
ity is required to be x 2 < 1-0. The transverse 
momentum vector of the B s candidate must 
be close to the displacement of the secondary 
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vertex from the primary vertex: the cosine of 
the opening angle a between the two vectors 
must fulfill cos(a) > 0.995. The significance 
of the B s candidate flight length l xy in the 
transverse plane is defined as l xy /cr xy (illus- 
trated in Fig. 2), where a xy is the error on 
the flight length with mean (<J xy ) = 120 y«m; 
we require l xy /a xy > 18.0. 
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Fig. 2. Decay length significance in the transverse 
plane for signal and background MC events. Both 
histograms are normalized to unity. 

In high-pj_ gluon-splitting events the bb 
quark pair moves closely together due to 
their boost, and the two decay vertices of the 
resulting 6-hadrons cannot be well separated 
in all cases. However, because of color re- 
connection, the hadronic activity around the 
dimuon direction is enhanced compared to 
the signal decay (where only one colorless B s 
meson decays). This is exploited in isolation 
requirements: The isolation I is determined 
from the dimuon transverse momentum and 
charged tracks with p± > 0.9 GeV in a cone 
with half-radius r — 1.0 around the dimuon 
direction as I = p ± ^/(p x ^ + Etrk bx|); 
we require I > 0.85. 

Figure 3 illustrates the mass resolution 
obtained on the signal MC event sample. 
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Fig. 3. Dimuon mass distribution in signal MC 
events. The curve is a fit of two Gaussians, the 
displayed parameters indicate the average mean and 
sigma. The histogram is normalized to unity. 

The distribution is fit with two Gaussians, 
the quoted width a = 36.0 MeV is deter- 
mined according to 

o N 2 a 2 + N 2 a 2 
~ N 2 + TV 2 

where a n = 32.1 MeV (a w = 60.2 MeV) and 
N n = 0.08 (N w = 0.03) are the width and 
normalization of the narrow (wide) Gaussian, 
respectively. 

Given the limited statistics of the back- 
ground sample, no event remains after the 
application of all selection requirements. 
However, the absence of correlation to the 
other selection requirements allows a factor- 
ization of the isolation and % 2 requirements 
from the other cuts in the determination of 
the total background rejection factor. The 
dominant sources of uncertainty on the signal 
(±25%) and background (±1™%) y ield are 
the statistical component of the background 
sample, the impact of the misalignment on 
the transverse flight length significance, and 
the assumption of factorizing cuts. 

The total selection efficiency for signal 
events is e s = 0.019 ± 0.002 and the back- 



ground reduction factor is £b — 2.74 x 10~ 7 . 
With this selection, the first 10.0 fb^ 1 of in- 
tegrated luminosity will yield ns — 6.1 ± 0.6 
signal events and ub = 13.8ti3.g background 
events in a mass window of tob, ±0.1 GeV. 
With this background estimate, the upper 
limit on the branching fraction is B(B® —> 
fi+fj,-) < 1.4 x 10~ 8 at the 90% C.L. 

6. Conclusions 

This study is limited by the size of the back- 
ground MC sample. In the future, it will 
include larger background samples and a de- 
tailed simulation of rare 6-hadron decays. 
The search for B® — > promises an 

interesting start-up analysis with the pos- 
sibility of setting tight constraints on the 
MSSM. With sufficient integrated luminos- 
ity, the precision measurement of the B® — > 
yit + /Lt _ branching fraction will set constraints 
on models of new physics. 
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